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Abstract. We analyze the left-right asymmetry in the semi-inclusive deep inelastic scattering (SIDIS) 
process without introducing any weighting functions. With the current theoretical understanding, we find 
that the Sivers effect plays a key role in our analysis. We use the latest parametrization of the Sivers and 
fragmentation functions to reanalyze the production process and find that the results are sensitive 
to the parametrization. We also extend our calculation on the K ± production, which can help us know 
more about the Sivers distribution of the sea quarks and the unfavored fragmentation processes. HERMES 
kinematics with a proton target, COMPASS kinematics with a proton, deuteron, and neutron target (the 
information on the neutron target can be effectively extracted from the 3 He target), and JLab kinematics 
(both 6 GeV and 12 GeV) with a proton and neutron target are considered in our paper. 



1 Introduction 



Single spin asymmetry (SSA) provides us with a power- 
ful instrument to investigate the internal structure of the 
' nucleon and its history can date back to the 1970s. In 
, the early 1990s, the E704 Collaboration reported the ob- 
' servation of a large left-right asymmetry in p^p — > ttX 
process pQ. This demonstrated that the transverse spin 
effect is significant even in high energies. In order to ex- 
plain the unexpected phenomenology, Sivers [2J first sug- 
gested a possible mechanism, the so-called "Sivers effect" 
today. But it was immediately criticized by Collins |3,4 , 
who proposed another mechanism known as the "Collins 
effect" now. Later, in Ref. [5], by considering the soft ini- 
tial state interactions, it was argued that the Sivers effect 
might be allowed. It was not until in 2002 that people 
began to realize that the final state interaction plays a 
crucial role in leading an SSA in the SIDIS process [6]. 
Then after considering the gauge links, it was found that 
the Sivers effect, or even the Sivers distribution could ex- 
ist [7] , and the Sivers distribution may have different signs 
in the SIDIS and the Drell-Yan processes. Despite of the 
early theoretical debates, some phenomenological analy- 
sis [HHH] attempted to explain the E704 data, and it was 
shown that the Sivers effect is important and other effects 
might be suppressed. But recently, an updated work [10] 
reported that the Collins effect is not strongly suppressed 
any more after a correction of a sign mistake. In all of 
these phenomenological works, TMD factorization were 
assumed, but we should be aware that the TMD factor- 
ization has not yet been proved for the pp — ► ttX process. 



Contrast to the complexity of the hadron-hadron col- 
lision process, where both the initial and final states are 
hadrons, the semi-inclusive deep inelastic scattering (SIDIS) 
process provides a cleaner and simpler playground for ex- 
ploring the nucleon structure. The azimuthal angle depen- 
dence of the cross section for this process has been system- 
atically studied in Ref. [TT] , where different structure func- 
tions were defined according to different azimuthal angle 
dependences. By multiplying different orthogonal weight- 
ing functions, we can isolate different structure functions 
from each other, and then extract the distribution or the 
fragmentation functions from relevant terms. For exam- 
ple, the Collins or the Sivers effect has a sin(<^ + <j/ s ) or 
sin(0^ — (f) e s ) Q modulation, respectively. Under the guid- 
ance, the HERMES [Hj and COMPASS [T3] Collabora- 
tions studied the sin(0^ + <f) l s ) and sin(</>^ — (f>g) asym- 
metries, and have confirmed the existence of the non-zero 
asymmetries. In the near future, JLab also plans a high 
precision measurement through the SIDIS process with 
a beam energy upgrading to 12 GeV. More such weight- 
ing SSAs will be studied in these experiments, and we 
hope these new observations will bring us more knowledge 
about the nucleon structure. 

If we turn back to the E704 observation, we find that 
the experiment just studied a simple un-weighted left- 
right asymmetry. The main reason is that in the inclu- 
sive hadron production, only one hadron is detected so 
that only one azimuthal angle can be defined. But in the 



1 (jihi^s are a l so written as cf>h,(f>s in some other literatures, 
but in this paper, we will write the explicit form with a super- 
script £ to address the lepton angle dependence. 
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SIDIS process, both the outgoing lepton and the produced 
hadron are measured, so we need to study a more com- 
plicated azimuthal dependence involving two azimuthal 
angles. However, we could still study the left-right asym- 
metry in the SIDIS process as the E704 experiment did. 
We suggest applying this method as an optional choice in 
analyzing the data, for it is a simple and basic quantity. 

In our previous paper [15] . we have studied this left- 
right asymmetry for the SIDIS process in the pion pro- 
duction. It was demonstrated that Sivers effect plays the 
most important role in producing a left-right asymmetry 
for a SIDIS process. In this paper, we will update the cal- 
culation with the new parameterizations of the DFs and 
FFs, and extend the calculation to the K ± production. 
As we know, the contribution from sea quarks, especially 
the ss quarks, might not be ignored in the kaon produc- 
tion. Also, we will extend our calculation to more kinemat- 
ics and targets for our prediction. HERMES kinematics 
with a proton target, COMPASS kinematics with a pro- 
ton, deuteron, and neutron target (extracted from the 3 He 
target), and JLab kinematics (both 6 GeV and 12 GeV) 
with a proton and neutron target are all considered in our 
paper. The main purpose of this paper is to reproduce a 
left-right asymmetry under the current theoretical frame- 
work, rather than to study the Sivers effect extensively. 
In our calculation, we will use the TMD distribution and 
fragmentation functions, as they may provide us a more 
vivid 3-dimensional picture of a nucleoli. The proof of the 
TMD factorization for a SIDIS process can be found in 
Ref. [H]. We will present our calculation up to a leading 
order approximation. 



2 Definition of the asymmetry and a 
theoretical description 

Obviously, the asymmetry is a function of space coordi- 
nates, so it depends on the choice of the coordinate sys- 
tem. For the experiments, the most convenient choice is 
to choose the direction of the beam and the target polar- 
ization as the spin plane. We call it as the £p frame, i.e. 
the laboratory frame, in which we can identify left or right 
according to the spin plane. We might perform the mea- 
surement in a space region left to the spin plane, then in 
its mirror region on the right, and by their differences, we 
can define a left-right asymmetry which is what the E704 
experiment did. We have noticed that changing the de- 
tected region from left ro right is equivalence to reversing 
the target polarization. We can express the asymmetry as 



.4 



1 Nj^Ps) - N{^js + tt) 
St N(ips) + N(tps + t) 



1 da^ — da- 1 
St da" + da^ 



where St is the transverse polarization of the target, and 
tps is the azimuthal angle of the spin vector. The minus 
sign in front of the expression is due to the fact that the 
detection occurred to the right of the beam, the same as 
that in the E704 analysis, where no weighting function is 
multiplied. 



Before our calculation, we first give an explanation 
to our kinematics. For a theoretical description, the £p 
frame is not always convenient, since we usually regard 
the SIDIS process as a virtual Compton scattering. So it 
is convenient to choose the "f*p frame, in which the z axis 
is defined along the direction of the exchanged virtual pho- 
ton, and the spin plane is defined by the virtual photon 
and the spin vector. Thus we have two reference frames, 
the £p and the j*p frames. In the Ip frame, we can de- 
fine the transverse spin vector St , the azimuthal angle for 
the spin vector and the produced hadron as ips and iphi 
respectively. All the kinematics can be manipulated eas- 
ily in this frame, but for the theoretical description, the 
7*p frame might be more convenient. In the j*p frame, 
we could define <fi e , (f>s and <j>h as the azimuthal angles for 
the lepton plane, the spin plain and the produced hadron 
plain, with respect to the horizontal plain in the labora- 
tory. Then we will define 4> e h = (f>h — <fi e , 4> l s = 4>s — <fi e , 
and these two angles are consistent with the Trento con- 
vention [16J 2 ! . 

The explicit expression for a SIDIS process can be 
found in Ref. [Hl[l~7] . where all the coordinate variables 
are defined in the "f*p frame. The connection between the 
two frames is via a rotation by a angle 9, due to which, 
a transverse spin vector in the laboratory frame has a 
longitudinal projection along the virtual photon [T51ll9| . 
Generally, 9 is very small and we will make further discus- 
sion later. By taking into account this, the cross section 
can be written as |19j : 

da 



dxdydtfig dzdtfi e h dP 2 ± 

^,2 



2sa;(l-e) 1 - sin 2 9 sin 2 0| 
St cos 9 



1 - sin 2 9 sin 2 i 



: sin 



St cos 9 



1 — sin 2 9 sin 2 i 



: sin 



'h - vs 



'h T VS 



M p 
h ■ fe_i 



JitDi 



+ other terms 
= c&ttju + dasw + dacoi + ■ ■ ■ , 
where we use a compact notation: 



(2) 



T[ufD\ = J> 2 f d 2 P± d 2 k ± 5 2 ( P± -k±- P h Jz) 

u( Px ,k ± )f a (x,p 2 x )D a (z,z 2 ki), (3) 

and 

i - y - \y 2 i 2 



i-y + i)y 2 + jy 2 j 2 '' 



h=P h± /\P h± \. (4) 



2 These two azimuthal angles are denoted as (f>h and 0s in 
the Ref. [16], but the notations are for other uses in our paper. 
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Table 1. An estimation on sin 9 





HERMES 


JLab 






6 GeV 


12 GeV 


(x) 


0.09 


0.23 


0.23 


(y) 
(Q 2 ) 


0.54 


0.6 


0.57 


2.41 GeV 2 


1.8 GeV 2 


2.5 GeV 2 


(sin 9) 


0.073 


0.19 


0.17 



First, we make a first approximation that 9 is small 
(We will give a detailed discussion later), thus the tp and 
7*p frames are of no difference. Now under this approxi- 
mation, we can use the kinematics defined in the j*p frame 
instead to analyze the asymmetry defined in Eq. [T] 



Aut(%,V, z,Ph±) 



1 /- 



b l h (das™ + da Co \ + •••), 



/ d(/) e s d4> l h dcruu 



f5) 



Next, we change the integral measure from d<p s d(f> e h to 
d(f> e d<f)h (The jacobian \ J\ = 1), and perform the integral 
over <j> e . We notice that sin0^ explicitly depends on <fr 
and all the convolution integrals T are independent of <j> . 
For the azimuthal angle dependence, we find that all the 
factors are oscillation functions of <// except sin(0/j — </>s), 
which is (/^-independent. So after integrating over <jr, we 
find that Sivers effect is o(l), but other terms such as 
the Collins effect are o(sin 2 9) (See detailed discussion in 
Ref. [T5]). So Sivers effect is dominant and other effects 
are suppressed in our analysis. We could also understand 
this result as the following. If we ignore the angle 9, we 
could set the 7*p frame equal to the Ip frame in the labo- 
ratory. In the E704 method, only one azimuthal angle was 
involved in fact, i.e. 4>h in our notation (not and the 
lepton angle <j/ in not included in the analysis. For a SIDIS 
process, only the Sivers effect is independent of the lepton 
plane, so it is not strange that the Sivers effect plays the 
most important role in leading an left-right asymmetry. 

We could make an estimation on the effect resulted 
from the angle 9. This angle can be calculated from the 
kinematical variables [19] 



smfl = 7-1 



i - y - \v 2 i 2 



1 + 7 2 



7 = 2xM p /Q. (6) 



We replace each variable by its average value to estimate 
the mean value of sin 9 for the HERMES and JLab exper- 
iments. 

The estimated result is shown in Table [TJ and we find 
that for most instance, the direction of the virtual photon 
is very close to the direction of the incident beam. There- 
fore, for convenience, the Collins effect which is not known 
so clearly yet is not considered in our analysis. 



3 Parametrization for distribution and 
fragmentation functions 

As a preparation for our calculation, we will present the 
parametrization for the distribution and fragmentation 
functions we will use in this section. 

For the Sivers functions, there are already some model 
calculations [50], but we would use a phenomenological 
parameterization for the Sivers functions. We should be 
cautious that a universal transverse momentum depen- 
dent Sivers distribution for different processes does not 
exist [21]. Fortunately, we will calculate for the SIDIS 
process, and the parametrization of the Sivers function 
is also from the SIDIS data. The Sivers effect has al- 
ready been studied by HERMES and COMPASS Collab- 
orations, and extractions on the Sivers functions for the 
u and d quarks were already obtained 8,22,23,24 . But 
all these results were under low statistics and assumed 
the existence of a symmetric and negligibly small Sivers 
sea. Recently, the HERMES Collaboration has provided 
much higher statistic data on the Ay^' 1 "^^ azimuthal 
asymmetry [55]. Besides the charged pion production, neu- 
tral pion and charged kaon azimuthal asymmetries were 
also analyzed. Also, the COMPASS Collaboration sep- 
arated the charged pion and charged kaon asymmetries 
from the charged hadron production measurement [? I26| . 
These SIDIS experimental data on the Sivers asymmetries 
for the pion and kaon production give us an opportunity 
to study the sea-quark Sivers functions for the u, d, s, and 
s quarks. With these data, in Ref. [57], the extraction of 
these functions was improved and the first estimates of 
the sea-quark Sivers functions were presented. The Sivers 
function is parameterized in the form 

= -^M q {x)f q {x)g{pl)h{pl), (7) 



Af q (x) = N q x^l-xf q ^ + ^ . 

a q p q 



(8) 



2 e-Pi/<fi> 2 
3(p_l) = / 2\ » h (P±) 



E± p -p1/(m 12 ) 



<2e- — e 



(9) 



All the parameters can be found in Ref. [57]. In the 
above parametrization, f(x) is the unpolarized parton dis- 
tribution functions, and we adopt the CTEQ6L parametriza- 
tion [55J as an input. We plot the Sivers functions for dif- 
ferent quark flavors in Fig. [1] 

For the fragmentation functions, all the former analysis 
of the fragmentation functions were based exclusively on 
the single-inclusive e + e _ annihilation (SIA) data and have 
been chosen the most simple functional form NiZ ai (\ — 
z)P* to parametrize the . But in these experiments, 
information on the quark and anti-quark fragmentation is 
always combined, for it always refers to the charge sum for 
certain hadron species, e.g. 7r + +7r~. In order to distinguish 
"valence" from "sea" fragmentation, some assumptions 
were proposed, e.g., in Ref. [29], D\ /D* = (1 — z) was 
assumed. In the last few years several one-particle inclu- 
sive measurements coming from both the proton-proton 
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collisions and the deep-inelastic lepton-nucleon scattering 
gave an opportunity to weigh each quark contribution in 
the hadronization process. In Ref. [3D], a global analysis 
was taken for the first time to analyze the individual frag- 
mentation functions for all flavors as well as gluons. A 
more flexible input is used 

H A^(l~z)ft[l + 7t (l-z)^] 

1 [Z ' Mo; B[2 + ai , 0i + 1] + 7i B[2 + a h A + Si + 1] ' 

(10) 

where B[a, b] is the Beta- function and iVj is normalized 
to represent the contribution of Df to the sum rule. For 
the fragmentation to 7r + , the isospin symmetry for the sea 
fragmentation functions is imposed, i.e., 

Df=Df. (11) 

But slightly different normalization in the q + q sum is 
allowed: 

D dU=NDZ + +u . (12) 
For the strange quarks it is assumed that 

D? =Df =N'Df. (13) 

For the charged kaons, it is assumed that the unfavored 
fragmentation functions are the same, 

=Df =DT =Df. (14) 

^From these relations, the unfavored fragmentation func- 
tions can be distinguished form the favored ones. Detailed 
parametrization for the integrated fragmentation function 
D(z) can be found in Ref. [30] . We present the numeri- 
cal results for the charged pion and kaon fragmentation 
functions in Fig. [5] and Fig. [3] We notice that in this 
parametrization, the s — > K + process is the most favored 
one, which means that sea quarks, especially s quark, 
might contribute significantly, although sea quark distri- 
butions are small compared with the valence quarks. So 




Fig. 3. Fragmentation functions for K + according to Ref. [30] 

measurements on the kaon production may help us to 
know more about the s(s) distribution. In our calcula- 
tion, we need the TMD fragmentation function, and we 
adopt a Gaussian assumption 

A^^^^ ^-y , as) 

with R 2 = 0.2 GeV 2 suggested in Ref. [31]. 



4 Numerical calculations 

The kinematical cuts used in the calculation are shown 
in Table O For the HERMES experiment, only the pro- 
ton target is calculated. For the Compass experiment, the 
proton, neutron and deuteron targets are all considered, 
while for the JLab experiment, the proton and neutron 
targets are assumed. In fact there is no free neutron tar- 
get, and in experiments the polarized 3 He is used. The 
effective asymmetry on a free neutron can be extracted 
from a 3 He target asymmetry. Detailed discussions can be 
found in some theoretical works [35] and a JLab's pro- 
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Table 2. Kinematics 





HERMES 


COMPASS 


JLabl 


JLab2 








proton 


neutron 


proton 


neutron 


Pbcam/GeV 

Q 2 /GeV 2 
W 2 /GeV 2 


27.6 


160 


6 


6 


12 


12 


> 1 


> 1 


> 1 


1.3 ~ 3.1 


> 1 


> 1 


> 10 


> 25 


> 4 


5.4 ~ 9.3 


> 4 


> 2.3 


X 


0.023 ~ 0.4 




0.1 ~ 0.6 


0.13 ~ 0.4 


0.05 ~ 0.7 


0.05 ~ 0.55 


y 


0.1 ~ 0.85 


0.1 ~ 0.9 


0.4 ~ 0.85 


0.68 ~ 0.86 


0.2 ~ 0.85 


0.34 ~ 0.9 


z 


0.2 ~ 0.7 


0.2 ~ 1 


0.4 ~ 0.7 


0.46 ~ 0.59 


0.4 ~ 0.7 


0.3 ~ 0.7 



„?0.1 



0.05 


-0.05 



Proton 






\ 

\ 

\ 



0.2 0.4 0.6 0.2 0.4 0.6 O.t 



Fig. 4. The x and 2-dependence of the left-right asymmetry 
for ir ± production on HERMES kinematics. Solid lines for 7r + 
and dashed lines for tt~ . Thick curves are our results and thin 
curves are results from Ref. [15] 



^0.05 



Proton 












Neutron 




\ \ 




Deuteron 









Proton 






s 


Neutron 









0.6 0.2 0.4 0.6 0.8 

Z 



Fig. 6. Similar as Fig. [4] but at JLab kinematics with a beam 



. GeV. 


Proton 








Neutron 













0.2 0.4 0.6 0.2 OA 0.6 

X Z 



W' 10- 2 10- 1 0.2 0.4 0.6 0.8 

x z 
Fig. 5. Similar as Fig-H but at COMPASS kinematics. 



Fig. 7. Similar as Fig. [4] but at JLab kinematics with a beam 
energy of 12 GeV. 



posal [33J . We will investigate the x and z dependence^ of 
the asymmetries. 

First, we present our results on -n^ production at dif- 
ferent kinematics. Fig. 5] - Fig. [7] show the results of the 

3 The E704 experiment only showed the dependence on if, 
i.e. approximate z here. 



left-right asymmetry. In Fig. [4] we make a comparison 
with the results already obtained in Ref. [TS], and we 
find that the two results are a little different. For the x- 
dependence of the asymmetry, our new results are sup- 
pressed when x increases. We can find the reason from 
the parametrization for Sivers functions as Fig. [1] shows. 
The new parametrization shows that at large x region 
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0.1 







-0.1 



-0.2 



Proton 








Neutron 









X Z 

Fig. 11. Similar as Fig- [HI but at JLab kinematics with a beam 
energy of 12 GeV. 

(x > 0.3), the u quark distribution falls down much faster 
than a previous parametrization, while the d quark distri- 
bution does not change so much. Notice that the u quark 
distribution gives a positive contribution, then we could 
understand our new results. For the z-dependence, the 
difference results not only from the distributions, but also 
from different parametrizations of the fragmentation func- 
tions. So we can say that our results are sensitive to the 
parametrization. 

Next we will extend our calculation to the K pro- 
duction, and the results are shown in Fig. [8]- Fig. QT] The 
K + production is quite similar to that of 7r + , but for the 
K~ production, we should be cautious. K~ is made up of 
a u and an s quark. So the sea quark contribution (from u 
and s) might be enhanced due to the favored fragmenta- 
tion process, while the valence quark contribution would 
be suppressed. Thus, the K~ production is a good way to 
study the sea quark distributions and the unfavored frag- 
mentation processes. Notice that the Sivers distributions 
of the u and s quarks given in Ref. [27] are so small with a 
large uncertainty that even their signs are not determined 
within the error, so our prediction must cover a large un- 
certainty area. Nevertheless, we hope that our results will 
be helpful to the future experiments and we expect that 
further experiments with a high precision will clarify the 
detail. 



5 Conclusion 



6 



Proton 









0.2 0.4 0.6 0.2 0.4 0.6 0.8 

X Z 



Fig. 8. The x and ^-dependence of the left-right asymmetry 
for K production on HERMES kinematics. Solid lines for K + 
and dashed lines for K~ . 
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-0.1 



Proton 








Neutron 








Deuteron 








\ 





10~ 3 10- z 10- 1 0.2 0.4 0.6 0.8 

x z 



Fig. 9. Similar as Fig. [8] but at COMPASS kinematics. 



o 



Proton 








Neutron 









X Z 

Fig. 10. Similar as Fig.JS] but at JLab kinematics with a beam 
energy of 6 GeV. 



Following the E704 analysis, we reanalyzed the left-right 
asymmetry in the SIDIS process with the new Sivers func- 
tions and fragmentation functions. In this paper, we con- 
sidered all the flavor contributions, including the sea quarks. 
We extended our analysis to the production process, 
and meanwhile to various kinematics with different tar- 
gets. We found that our results are sensitive to the parametriza- 
tion form of the distribution and fragmentation functions, 
so we consider it necessary to perform higher precision 
measurements to constrain the parametrization. 
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Our prescription originated from the E704 experiment 
is an optional and simple way analyzing the data. In this 
prescription, no weighting functions are multiplied, al- 
though it might not give any more information. We sug- 
gest that relevant experiment collaborations could present 
their data in this new way as an optional choice for further 
theoretical studies. 
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